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Abstract 

Using the ZEUS detector at the HERA electron-proton collider, we have searched 
for lepton flavor violation in ep collisions at a center-of-mass energy of 300 

GeV. Events of the type e + p — > £ + X with a final-state lepton of high transverse 
momentum, t = /i or r, were sought. No evidence was found for lepton flavor 
violation in the combined 1993 and 1994 data samples, for which the integrated 
luminosities were 0.84 pb~ l for e~p collisions and 2.94 pb~ l for e + p collisions. Limits 
on coupling vs. mass are provided for leptoquarks and i?-parity violating squarks. 
For flavor violating couplings of electromagnetic strength, we set 95% confidence 
level lower limits on leptoquark masses between 207 GeV and 272 GeV, depending on 
the leptoquark species and final-state lepton. For leptoquark masses larger than 300 
GeV, limits on flavor-changing couplings are determined, many of which supersede 
prior limits from rare decay processes. 
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1 Introduction 

Lepton flavor is conserved in all interactions of the Standard Model (SM); discovery of 
lepton-flavor violation (LFV) in any form would be evidence for physics beyond our 
principal particle physics paradigm. Many searches for specific reactions which violate 
lepton flavor have been performed. The most sensitive include searches for \i + N — > 
e + N using very low-energy muonsfll]], for the forbidden muon decay /i — > e-y 0, and 
for forbidden leptonic kaon decays The limits from these processes are sensitive to 
e fi flavor change, but not to e <-» r. Also, each of these processes involves specific 
quark flavors: in the first case, only first generation quarks participate; in the second case, 
for mechanisms which involve virtual quarks, the same quark flavor must couple to both 
e and /i; in the last case, strange quarks must be involved. Since lepton flavor change 
could involve the r-lepton or could be accompanied by quark flavor change, there may 
be LFV reactions which would be invisible to these very sensitive experiments. Hence, 



we have no a priori reason to assume that flavor violation, should it exist, will be visible 
through these specific reactions. Therefore, though less sensitive in an absolute sense, 
other manifestations of flavor violation, like forbidden leptonic decays of B- and D— 
mesons and of r-leptons, are being investigated ||] . 

We report here a search for LFV carried out by the ZEUS collaboration at the HERA ep 
collider where we have sought instances of the reaction 

e + p^i + X, (1) 

where I represents an isolated final-state \x or r with large transverse momentum and 
X represents the hadronic final state. Processes with such topologies can be found in 
ZEUS with good efficiency and with little background. It should be emphasized that any 
reaction of the type ([!]) in which a final-state high-energy /iorr replaces the incident 
electron^ would be direct evidence for physics beyond the Standard Model, independent 
of the underlying mechanism. Furthermore, this reaction should occur at some level for 
a wide range of possible LFV mechanisms. LFV mechanisms which also involve a quark 
flavor change, or which are stronger for heavier quarks || may be seen more readily at 
HERA, where the sensitivity is largely independent of quark flavor]^, than in low-energy 
experiments. 

The lepton-flavor violating reaction ([I]) could occur via s-, u—, or t-channel exchanges 
as shown in figure [l|. For the s— and ^-channel processes the exchanged particle has the 
quantum numbers of a leptoquark (or an i?-parity violating squark). The cross sections 
depend on the leptoquark species and mass, and on the couplings, X eqi and \e q2 shown 
in figure [I]. For the case of t-channel exchange, the process would be mediated by a 
flavor-changing neutral boson. 

For definiteness, we will describe reaction (H) with leptoquarks as the carrier of the LFV 
force, treating separately the cases of direct leptoquark production and the virtual effects 
of leptoquarks with masses above 300 GeV. The similarity of production formulae between 
i?-parity-violating squarks and certain leptoquarks permits us to relate the couplings 
implied by the two mechanisms for a specified cross section. Results on flavor violation 
induced by leptogluons or by flavor-changing neutral bosons as well as details of the 
technique used in this analysis are also available 0. The HI collaboration [|7| has also 
searched for direct production of leptoquarks with flavor- violating couplings using similar 
methods. 

This analysis is based on an integrated luminosity of 0.84 pb^ 1 (2.94 pb^ 1 ) of e~p (e + p) 
data taken during the 1993 and 1994 running periods. Beam energies at HERA were 26.7 
GeV (27.5 GeV) for the electron beam in 1993 (1994) and 820 GeV for the proton beam. 
The resulting center-of-mass energy of 296 GeV (300 GeV) is an order of magnitude 
higher than for fixed-target lepton-nucleon scattering experiments. 



1 In the following, "electron" is generically used to denote both electrons and positrons. 

2 Though the threshold for top (t) quark production is below the HERA center-of-mass energy, with 
present luminosities t production would only be observable if the couplings were very large. Therefore, 
we choose not to report on LFV couplings involving top in this paper. 
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2 Scenarios of Lepton Flavor Violation 



2.1 Leptoquarks 

A leptoquark (LQ) is a hypothetical color triplet boson with fractional electric charge, 
and non-zero lepton and baryon numbers. Such particles are often invoked in extensions 
of the SM, e.g. in grand unified theories and technicolor models f|. It is possible, indeed 
desirable in some models, that a LQ couple to multiple lepton and/or quark flavors, 
thereby providing a mechanism for flavor violation. 

The simplest models involving a flavor-violating leptoquark would be characterized by 
three parameters: the leptoquark mass, and the coupling at each lepton-quark-leptoquark 
vertex. In order to avoid models which would involve additional parameters, we have 
assumed the following four points: 



1) the LQ has SU(3)c X SU(2)i x U(l)y invariant couplings, 

2) the LQ has either left- or right-handed couplings, but not both (i.e. XlXr = 0), 

3) the members of each weak-isospin multiplet are degenerate in mass, 

4) one LQ species dominates the production process. 



There are fourteen species of leptoquarks which satisfy these conditions |§ . For fermion 
number F = L + 3B (L and B denote lepton and baryon number) equal to zero, the 
species are denoted [|] Sfy 2 , 5*f /2 , S^ /2 , V L , V R , V R , and Vf. For F = 2, they are Sfr, 
Sq, Sq, Si, Vy 2 , Vy 2 , and Vy 2 . Here S and V indicate scalar and vector leptoquarks 
respectively, which couple to left- (L) or right-handed (R) leptons as indicated by the 
superscript. The subscript gives the weak isospin of the LQ0. In s-channel reactions, 
F = LQ cross sections are higher in e + p collisions, where they are produced via e + q 
fusion, than in e~p collisions where e~q fusion occurs. The reverse is true for an F = 2 
leptoquark. 

A LQ scenario is defined by the leptoquark species, by the generations of the quarks 
which couple to the electron and to the final-state lepton, and by the final-state lepton 
flavor. Hence there are 14x3x3x2 = 252 different LQ scenarios, each characterized by 
two dimensionless couplings, X eqi and A^ 92 , defined in figure p], which could induce flavor 
violation. Such LQs would also mediate flavor-conserving interactions with a final-state 
e or u e , which are not considered in this paper. 



As an illustration, we show in figure |2] the present limits |T0[ on \ eqi \£ q2 versus LQ 
mass, M LQ , for reactions which could proceed through the left-handed scalar isosinglet 
LQ, Sq. Note that each of the limits assumes that specific quark flavors couple to e 
and [i. For example, the most sensitive limit, from \xN — > eN, applies only for first- 
generation quarks in both initial and final states. Also in this figure are the results of 



direct searches for leptoquark pair production. Searches in e + e collisions at LEP[1I 



3 The tilde differentiates LQ species which differ only in that one species couples to u-type quarks and 
the other to d-type quarks. See M for details. 
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exclude scalar leptoquarks lighter than about 45 GeV which couple to e, /x, r, or any 
neutrino. For leptoquarks with couplings of electromagnetic strength, masses below 73 
GeV are excluded. While the LEP experiments did not search for flavor violation, their 
non-observation of eeqq, fijlqq, TTqq, or vVqq final states kinematically consistent with 
leptoquark pair production imply flavor-violating leptoquark mass limits which are weaker 
by at most a few GeV. At the Tevatron []T2J , searches in pp collisions have excluded scalar 
leptoquarks lighter than 131 GeV (96 GeV) for an assumed branching fraction to q\i of 
100% (50%). These limits on M LQ are independent of the LQ couplings in most models. 

The LQ-induced cross sections for reaction ([[J), given in detail in the appendix, depend 
on the initial quark density, the couplings, and the species of LQ involved in the reaction, 
as well as on the kinematic event variables x (the Bjorken scaling variable) and y (the 
inelasticity). Here x is defined as x = —q 2 /(2q ■ P) and y as y = (q ■ P)/(k ■ P), where 
k, k', and P are the four-momenta of the initial-state electron, the final-state lepton 
and the proton respectively, and q = k — kl '. The square of the center-of-mass energies 
of the electron-proton and the electron-quark systems are given by s = (k + P) 2 and 
s = xs respectively. The remaining Mandelstam variables are given by t = —sxy and 
u = — sx(l — y). 

For a given coupling, the cross section is largest when the LQ mass, M LQ , is less than 
yfs. In this case, the LQ is produced in the s-channel, as indicated in figure [l|a. Such a 
leptoquark will appear as a narrow resonance in the ^-distribution peaked at Xq = M 2 Q /s. 
In the narrow-width approximation described in the appendix, the cross section for this 
process using unpolarized beams can be written 

a eqi ^ £q2 = ^-X 2 eqi B eq2qi (x ,Ml Q ) dyj{y), (2) 



where f(y) 




scalar LQ 
y) 2 vector LQ, 



and qi(x, M 2 Q ) is the quark density in the proton for the initial-state quark (or antiquark) 
flavor qx, X eqi is the coupling at the LQ production vertex and B^ q2 is the branching fraction 
of the LQ to lepton I and quark flavor q^. In this process, the final state lepton will have 
a transverse momentum (P/) of order M LQ /2. 



For resonant s-channel production, the cross section for flavor-violating events is propor- 
1 n We will set limits on this quantity as a function of M LQ . 



tional to \ eq J5t q2 



For the case M LQ 3> y/s, either or both s- and it-channel contributions may be important. 
The corresponding cross sections can be written as 

n 2 



eqi—>lq2 



eq2^lqi 




where f(y) 



1(1 
2(1 

2 



y? 
yf 



dx dyxq 1 (x,s)f(y), 
dx dyxq 2 (x, -u)f(y), 



s-channel scalar LQ 
li-channel scalar LQ 
s-channel vector LQ 
w-channel vector LQ, 



(3) 
(4) 



o 



and the indices q\ and q2 specify the quark flavors which couple to the electron and the 
final-state lepton, respectively. Here the final-state lepton again will have large transverse 
momentum with Pf « vI/2. 

Notice that in the high-mass case, all information about the leptoquark mass and cou- 
plings is contained in the quantity X eqi \e q2 / M% Q , which is the quantity on which we set 
limits. As might be anticipated, other LFV processes mediated by leptoquarks, such as 
flavor violating meson decays, are sensitive to exactly this quantity. Hence, our results 
may be compared directly with prior LFV searches. This is done in section |5]. 



2.2 R— Parity Violating Squarks 

Squarks (q) are the hypothesized supersymmetric partners of quarks. In supersymmetry 
(SUSY), _R-parity is defined as Rp = (— 1} 3B + L + 2S where B, L, and S denote baryon 
and lepton numbers and spin respectively. This implies that Rp — +1 for SM particles 
and Rp = — 1 for SUSY particles. If i?-parity were conserved, SUSY particles would be 
produced in pairs and ultimately decay into the lightest supersymmetric particle (LSP), 
which would be stable and neutral. We refer here to this LSP as the photino (7). In 
a model with i?-parity violation, denoted Rp, single SUSY particle production would 
occur and the LSP would decay into SM particles. Of particular interest for ep collisions 
are i?-parity violating superpotential terms of the form |13| X^L^Q^D R . Here L L , Q L , 



and Dp denote left-handed lepton and quark doublets and the right handed <i-quark 
singlet chiral superfields respectively, and the indices i, j, and k label their respective 
generations. Expanded into four-component Dirac notation, the corresponding terms of 
the Lagrangian are 

£ = K jk + d L d k R ui + (4)*(4) c 4 - ~ %?rA ~ (4)*(el)X] + h. c. 

(5) 

For i — 1, the last two terms will result in u and d production in ep collisions. Identical 
terms are found in the Lagrangians for the scalar leptoquarks S1/2 and So, respectively 
0- 



Lepton-flavor violating ep interactions would occur in a model with two non-zero couplings 
A^-fc which involve different lepton generations. For example, the process ed — > tt J — > ~pd k 
shown in figure ^|a involves the couplings and \' 2 j k . Similarly, non-zero values for 
A' llfe and \' 3 j k would lead to the reaction eu — > d k — >• ru^ shown in figure [|b. Down-type 
squarks have the additional decay d k — > v l d? , a mode unavailable to up-type squarks. 

The difference between mechanisms involving i?-parity violating squarks and leptoquarks 
is that the squarks may have additional i?-parity conserving decay modes with final-state 
neutralinos, such as q — ► q^ (shown in figure |3]c) or with final-state charginos, as in 
u — > dx + ■ The branching ratios B q ^ for the i?p-conserving decay q —> qj and B'^ k for 
any Rp decay mode are related |L|] by 



R' R 

where A^- fc is the Rp coupling at the decay vertex, a is the electromagnetic coupling^, e q 
4 We evaluate a at the scale Mz (a=l/128) because s is of order Mz at HERA. 
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is the squark charge in units of the electron charge and the photino and squark masses 
are and rriq, respectively. 

Coupling limits for LFV decays of an Sq leptoquark can be interpreted as d k coupling 
limits through the correspondence X eqi yj Bt q2 = X' nks jB~l~^, where i and j are the genera- 
tions of the LQ decay products I and q 2 . Similarly, coupling limits on the Sy 2 LQ can be 
converted to limits on couplings to vP via X eqi \jB^ q2 = X^^^J B[- k} where i and k are the 
generations of £ and q^. 

If the stop (t) (15| is lighter than the top quark, then the i?p-conserving decay t — > £7 
(figure |3|c) will not exist. In the case of t, the correspondence with the coupling limit 
on Sy 2 is given by X e d^J Bg q2 = cos t X' 131 ^J B' i3k where 6 t is the mixing angle between the 
SUSY partners of the left- and right-handed top quarks. Over a broad range of possible 
stop masses, it is expected that cos 2 6 t ~ 0.5 fl5fl . 



3 The ZEUS Detector and Event Simulation 



The main components of the ZEUS detector 16 used for this analysis were the uranium- 



scintillator calorimeter (CAL) |l7j and the central tracking detector (CTD) fl 

The CAL, which covers polar angles^ between 2.2° and 176.5°, is divided into forward 
(FCAL), barrel (BCAL), and rear (RCAL) parts. Each part is further subdivided into tow- 
ers which are longitudinally segmented into electromagnetic (EMC) and hadronic (HAC) 
sections. In depth, the EMC is one interaction length; the HAC sections vary from six 



to three interaction lengths, depending on polar angle. Under test beam conditions JT7 ] 



the calorimeter has an energy resolution of <7E(GeV) = 0.18 JE(GeV) for electrons and 



cr^GeV) = 0.35 yE(GeV) for hadrons. In this analysis, only cells with energies above 
noise suppression thresholds (60 MeV for EMC, 110 MeV for HAC) were used. 

A superconducting coil located inside the CAL provides a 1.43 Tesla magnetic field parallel 
to the beam axis in which the charged particle tracking system operated. The interaction 
vertex is reconstructed with a resolution of 4 mm (1 mm) along (transverse to) the beam 



direction. The muon detection system [19] was used to check the efficiencies and the 



background estimates for the primary muon identification, which used only the CAL and 
the CTD. The muon detectors are also divided into three sections covering the forward, 
barrel, and rear regions. In the barrel and rear sections, which were used for this analysis, 
the detectors consist of eight layers of limited streamer tubes, four layers on each side 
of the 80 cm thick magnetized iron yoke. Luminosity was measured p(| from the rate 
of bremsstrahlung events (ep epy) detected by a photon calorimeter (LUMI) located 
downstream of the main detector. The luminosity is known to 3% for the e~p data and 
to 2% for the e + p data. 

To evaluate detection efficiencies, we have simulated flavor- violating LQ processes using a 



5 The ZEUS coordinate system is right-handed with the Z axis pointing in the proton beam direction, 
hereafter referred to as forward, and the X axis horizontal, pointing toward the center of HERA. The 
polar angle 9 is taken with respect to the proton beam direction from the interaction point. 
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modified^ version of pythia [|2TJ and also with lqmgen which is based on the differential 
LQ cross-sections given in ||. The calculations included initial state bremsstrahlung. 



For background estimation, charged-current (CC) and neutral-current (NC) deep-inelastic 
scattering (DIS) events with electroweak radiative corrections were simulated using lepto 

via DJANGO 



22 interfaced to HERACLES 23 



The MRSA |25| parton density pa- 
rameterization was used. The hadronic final state was simulated using ARIADNE p6| and 



JETSET PT| . 

Photoproduction processes were simulated using herwig [P7[ , and photoproduction of cc 
and bb pairs by PYTHIA and AROMA f28fl . The processes 77 — > fi + [J>~ and 77 — > t + t~ 
were generated using zlpair ||29|1 . Finally, production of W bosons was simulated using 



EPVEC [|G 



All generated events were passed through a geant [[3TJ based detector simulation which 
tracked final state particles and their decay and interaction products through the entire 
detector. The simulated events were processed with the same analysis programs as the 
data. 



4 Trigger and Analysis 

The signature of LFV events (e + p — > £ + X) in this experiment is an isolated fi or r of 
high transverse momentum, P/ ~ v^S/2, balanced by a jet of hadrons. 



4.1 Search Strategy 

Our search strategy relies on the fact that the LFV signal events will almost always have 
a large net transverse momentum Pf al measured in the calorimeter. We reconstruct P£ al 
as ff = (Pf + P2)V2. Here Px = E . E . sin( 0.) cos( 0.) and Py = ^^sin^sin^) 
where the sums run over all calorimeter cells and E iy 8i, and 0j are the energy, polar angle 
and azimuthal angle of cell i, calculated using the reconstructed event vertex. We also 
reconstruct the azimuth of the missing transverse momentum, determined from cos <p miss = 
-Px/fr 1 and sin0_ s = -P Y / ' l . 

A high energy muon is a minimum ionizing particle, typically producing a measured energy 
of about 2 GeV in the calorimeter. If the much larger muon transverse momentum, P/*, 
is balanced by a jet of hadrons, then P^ al w P/\ Thus the signature for such an event 
would be a large P 4 ca/ and a high momentum track which points to an isolated calorimeter 
cluster with approximately 2 GeV of energy at an azimuthal angle 4> mi3S . 

A final-state r decays promptly to a small number of charged particles (1 or 3, 99.9% of 
the time), zero or more neutral hadrons, and at least one neutrino. Since the r mass is 
small (m T = 1.78 GeV) compared to its transverse momentum, the r decay products will 

6 The final-state electron and quark from these generators were replaced by the appropriate lepton 
(/1 or t), and quark species, before the simulation of parton showering and fragmentation. Both s- and 
it-channel exchange contributions were included. 
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be collimated in a cone of opening angle ~ 0.03 radians. For events in which the r decays 
via r — > \wv, the experimental signature will be similar to an event with a final state 
muon except with P/* < pf al . If the r decays via r — ► ei/77, the event will be characterized 
by large P^ al (due to the undetected neutrinos) and the presence of a high-transverse- 
momentum electron with azimuth <p mi8S . Finally, in the case of a hadronic r decay, we 
would again see a large P^ al due to the neutrino, and a compact hadronic cluster with 1 
or 3 tracks, also at azimuth 4> mi3S . 



4.2 Trigger 



Data were collected with a three-level trigger system |Tt|. Since the signature which we 
are seeking is one with missing transverse momentum measured in the calorimeter, our 
triggering scheme was largely calorimeter based. The first-level triggers used net trans- 
verse energy, missing transverse energy, as well as EMC energy sums in the calorimeter. 
The thresholds were well below the offline requirements. The second-level trigger rejected 
backgrounds (mostly p-gas interactions and cosmic rays) for which the calorimeter timing 
was inconsistent with an ep interaction. Events were accepted if P^ al exceeded 9 GeV 
and either a track was found in the CTD or at least 10 GeV was deposited in the FCAL. 
The latter alternative was intended to accept events with jets which are too forward for 
the tracks to be observed in the CTD. The third-level trigger applied stricter timing cuts 
and also pattern recognition algorithms to reject cosmic rays. 



4.3 Leptoquark Mass Reconstruction 

For M LQ < 300 GeV, the leptoquark is produced as an s-channel resonance and conse- 
quently, the invariant mass distribution of the q£ final state is sharply peaked at M LQ . 
When searching for a leptoquark of a given mass, the expected background can be reduced 
by requiring the reconstructed q£ mass to be consistent with M LQ . 

We reconstruct the leptoquark mass as follows using a simple ansatz based on three 
approximations: 1) the four-momentum of all final state muons and neutrinos can be 
represented by a single massless pseudoparticle; 2) the contribution of the proton remnant 
to the reconstructed mass can be ignored; and 3) no energy escapes through the rear beam 
hole. 

The four-momentum of the invisible pseudoparticle, P invis ; is related to the net four- 
momentum P = (E,Px,Py,Pz) measured in the calorimeter as P™ via = —Px, Pp- Vis = 
—Py, and E — Pz + E invia — p^ vis = 2E e where E e is the electron beam energy. The 
reconstructed leptoquark mass is given by (M[g) 2 = (P + p™"") 2 . 

We have applied this mass reconstruction to simulated LFV events and determined two 
functions, /^(M^q) and a r L g(M LQ ) which give the mean and the standard deviation of a 
Gaussian fit to the reconstructed mass distribution as a function of the true M LQ . Studies 
of simulated LQ signals indicate that the mass resolution improves from about 13% at 
M LQ = 100 GeV, to about 6% at M LQ = 250 GeV. 
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4.4 Event Selection 



The most important offline selection requires that P t ca exceed 12 GeV. The initial event 
selection is designed to accept all ep collisions which meet this condition, while efficiently 
rejecting the high-rate backgrounds from cosmic rays, proton-gas interactions, off-beam 
protons, and beam-halo muons. Triggers from these backgrounds usually do not have a 
reconstructed vertex. In cases where a spurious vertex is reconstructed, it typically is made 
from a small number of low-momentum spiraling tracks which do not intersect with the 
beam line. Unlike ep collisions, for which the distribution of Z vertex position is centered 
at Z = with an r. m. s. width of 12 cm, the spurious vertices have a Z distribution 
which is roughly uniform. In cases of protons colliding with residual gas in the beam pipe, 
or with the beam-pipe itself, the low-multiplicity spurious vertex is accompanied by a 
large number (10 to 100) of tracks which are not correlated with the vertex. Occasionally 
a cosmic ray or a beam-halo muon will coincide with an ep interaction which provides the 
reconstructed vertex. In these cases the vertex tracks are typically of quite low momentum 
(0(100 MeV)). 

In order to remove such backgrounds, we require that a vertex is reconstructed and that it 
lie within 50 cm of the nominal interaction point. We define N tr k to be the total number 
of reconstructed tracks, N goo d to be the number of tracks with transverse momentum 
P t > 300 MeV and a distance of closest approach to the beam-line of less than 1.5 cm, 
and N vtx to be the number of tracks forming the vertex. We require N gooc i > 1 and 
Ngood — 0.05iV fr .fc. In order to reject proton-induced background, for which the energy 
deposited in the calorimeter is concentrated at small polar angles, we remove events with 
P z /E > 0.8(0.94)[| if N trk - N vtx > 80 (20). In addition, we require the timing of each 
calorimeter cluster with energy above 2 GeV to be consistent with an ep interaction. To 
reduce the cosmic ray background, we apply an algorithm which rejects events in which 
the pattern of calorimeter energy deposits is consistent with a single penetrating particle 
traversing the detector. 

The 175 events which passed these cuts were visually examined and 29 events clearly 
initiated by cosmic rays, muons in the beam halo, or anomalous photomultiplier discharges 
were removed, leaving 146 ep collision events. These events were divided into two classes: 
those for which no isolated electron with energy E e > 10 GeV was found in the calorimeter 
(class and those for which such an electron was found (class e). The following selection 
cuts, which were developed in Monte Carlo studies, were applied to each sample in order 
to eliminate SM backgrounds. 

There were 114 events with no isolated electron. Six events were rejected because 
an electron of more than 5 GeV was observed in the luminosity electron calorimeter 
and they were thus recognized to be background due to photon-proton {^p) col- 
lisions. This left 108 a 1 events, which agrees well with the Monte Carlo estimates 
of 100 CC DIS events and 15 7p and 77 events. The Pjr al distribution of this $ 
sample is compared with the Monte Carlo prediction in figure |]a. The e' sample 
serves as the source of flavor violation candidates with a muon in the final state, as 

Here Pz and E are reconstructed from the calorimeter cells in a manner similar to the components 
Px and Py described above. 
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well as of candidates with final-state r's which subsequently decay via r — * [ivv or 
t — ► z/+hadrons. 

e: There were 32 events which contained an isolated electron. In order to reject NC 
DIS background, for which the electrons are concentrated at large polar angles, we 
required the electron polar angle to be less than 100°. After this cut, 12 events 
remained in the e sample, compatible with the Monte Carlo prediction of 14 NC 
DIS events. The f^ al distribution of these remaining events is compared with the 
Monte Carlo prediction in figure |]b. LFV candidates with a final-state r which 
decays via r — > evv (18% branching fraction) were sought in this sample. 

The final cuts rely on a clustering algorithm^] which assigns each calorimeter cell above 
noise threshold to one and only one cluster. Each cluster is characterized by its en- 
ergy, E clu , as well as the energy-weighted mean azimuth, (j) clu , and pseudorapidity, r) clu = 
— ln[tan(0 ciu /2)]. We expect the final-state lepton in a LFV event to produce a single 
isolated cluster. To decide if a cluster is isolated, we examine the set of all calorimeter 
cells which are within 0.8 units in f](f) of the cluster ({(4> c i u — (f> C eii) 2 + { r lciu — Vceii) 2 ] 1 ^ 2 < 0.8). 
A cluster is defined to be isolated if the summed energy of all calorimeter cells in this set 
which do not belong to the cluster is below 2 GeV. For each cluster, we also compute 4> clu , 
which is defined as the energy weighted mean azimuth of all cells in the entire calorime- 
ter, except for those assigned to the cluster. Note that for LFV events cf) clu differs only 
slightly from <f) miss which is computed using all calorimeter cells. A cluster is said to be 
ft al -aligned if it satisfies the inequality: cos(0 c;u — 4> clu ) < cos 170°. This ensures that the 
isolated cluster is opposite in azimuth to the rest of the energy in the calorimeter. 

To enter the final sample for \xq or rq final states, an event must satisfy the criteria of 
one of four selections, described below. 

H or r — > //: In a class event, there must exist an isolated /^-aligned cluster with 
energy 0.5 GeV < E clu < 6 GeV and at most 80% of its energy in the electromagnetic 
layer of the calorimeter. It must have exactly one matching trackf] and that track 
must have momentum exceeding 20 GeV. The efficiency^ to satisfy these cuts for 
scalar (vector) leptoquarks which decay to fiq decreases with LQ mass from 74% 
(78%) at M LQ = 80 GeV to 31% (50%) at M LQ = 260 GeV. The background 
estimate for the n selection was 0.1 events from the inelastic process 77 — ► fi + fi~. 
Zero events were observed in the data. 

To check the efficiency of the muon selection, we performed an independent event 
selection which did not use the CAL, but required a track in the barrel or rear muon 
chambers which was matched to a CTD track with a transverse momentum of at 
least 5 GeV. In order to select events with isolated muons, we rejected events which 
had an electron found in the calorimeter or had more than three tracks fitted to 

8 The clustering algorithm joins each cell with its highest energy neighbor, thus producing one cluster 
for each cell which has more energy than any of its neighbors. Two cells are defined as neighbors if they 
are in towers which share a face or an edge. Cells on the forward or rear edges of the BCAL are also 
neighbors with the FCAL or RCAL cells which are behind them, as viewed from the interaction point. 

9 A matching track is defined such that the distance of closest approach between the extrapolated 
track and the calorimeter cluster is less than 30 cm. 

10 All quoted efficiencies include the trigger efficiency. 
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the vertex. A total of 15 events, which contained 17 CTD-matched muon chamber 
tracks passed these cuts. This number agrees with the Monte Carlo estimate of 20 
events from 77 — ► All 17 tracks were matched to an isolated calorimeter 

cluster which passed the cuts described above (except for the /^-alignment). 

t — ► e: In a class e event, the isolated electron must be /^"'-aligned. The efficiency 
to satisfy these cuts for scalar (vector) leptoquarks with final-state rq, and the 
subsequent decay r — ► euu, rises with LQ mass from 23% (17%) at M LQ = 80 GeV 
to 75% (75%) at M LQ = 260 GeV. The background estimate for the r — > e selection 
was 0.2 events from NC DIS. Zero events were observed in the data. 

r — > hadrons: In a class $ event, there must exist an isolated /^'-aligned cluster with 
E clu > 10 GeV, which has either 1 or 3 matching tracks. At least one track must 
have a momentum exceeding 5 GeV. The efficiency to satisfy these cuts for scalar 
(vector) leptoquarks with final-state rq and hadronic r decay rises with leptoquark 
mass from 15% (12%) at M LQ = 80 GeV to 39% (47%) at M LQ = 260 GeV. We 
estimated a background of 0.4 events for this selection coming from CC DIS (0.2 
event), 77 — > (0.1 event), and cc production (0.1 event). We observed zero 

events. 

f t cal > 80 GeV: Leptoquarks with mass in the range 200 GeV< M LQ < 300 GeV would 
be strongly boosted in the forward direction so that the final state \x or r would 
often have polar angle less than 10°. In such cases, the final-state lepton would be 
outside the CTD acceptance and would consequently fail the track matching cuts. 
In order to maintain high efficiency at these masses, we accept any event from either 
class e or class for which f£ al > 80 GeV. For 240 GeV leptoquarks which decay 
to fiq, accepting events with f > ^ al > 80 GeV increases the overall acceptance from 
36% to 69% for scalars, and from 53% to 76% for vectors. For the f^ al > 80 GeV 
selection, we estimated a background of 1.0 events from CC DIS and we observed 
zero events. 

For the low-mass leptoquark search [M hQ < 300 GeV), one additional cut was applied, 
which, in contrast to all cuts described above, depends on M LQ , the mass of the LQ being 
searched for. The leptoquark mass was reconstructed using the method described in [4.3| 
and we required that the reconstructed mass must lie within 3a r L Q(M LQ ) of /i^(M £Q ). 



5 Results 

With no candidate events for LFV found in either the e~p or e + p data samples with 
integrated luminosities C e - = 0.84 pb -1 and £ e + = 2.94 pb~ x , we set upper limits on the 
couplings of the various LFV processes described in section |2|. The upper limit on the 
coupling A is obtained from the relation iV = J2i= e +,e- A^c^A), where e is the efficiency, 
er(A) is the cross section for a coupling A, and N is the Poisson 95% confidence level (CL) 
upper limit on the number of events. The signal efficiencies and background estimates 



were determined by Monte Carlo studies. We estimate the systematic uncertainties in 
the efficiencies to be 5%. Cross sections were calculated using the formulae given in 
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the appendix and the GRV-HO |33| parton density parameterization. Cross sections 
calculated using the MRSH [34] parameterization differ in magnitude by less than 12% 
for u, d, s, or c quarks in the initial state, and by less than 19% for initial-state b quarks. 



5.1 Low-Mass Leptoquark Limits (M^q < 300 GeV) 

In the case of low-mass leptoquarks, we calculate upper limits on X 2 B^ q2 using equation 
0. Since our limits are largely independent of the final state quark type (as long as the 
top quark is not involved), we show in figures |^ and |6] the upper limits on A egi J 1 B iq21 for 
scalar and vector leptoquarks where q\ is a first-generation quark and q 2 is the final-state 
quark of any generation (except top). The limits for t = \i and for i = r are shown 
separately as a function of LQ mass for the various scalar and vector LQ species. We note 
that for several LQ species, we probe coupling strengths as small as X 2 /An « 10~ 3 a for 
M LQ = 100 GeV and B lq2 = 0.5. 

In figure |7| we compare these limits on LFV with those from previous searches for two 
representative LQ species, Sq and V Q R . Assuming that Bg q2 = 0.5, we plot as a solid curve 
the upper limit on X eqi as a function of the LQ mass. Curves are shown for both \x (upper 
plots) and r (lower plots) final states. In contrast with many other limits on LFV, the 
coupling limits from this experiment apply to final-state quarks of any generation (except 
top). The various broken curves are low-energy limits quoted from reference For 
each of these curves, the pairs of numbers in parentheses denote the generations of quarks 
which couple to e and t. Coupling limits for Bi q ^ 0.5 can be obtained by multiplying the 

limit on X eqi plotted in figure [F] by ^J0.5/Bi q . We emphasize two important implications 
of figure [/J 



1. The ZEUS limits on ed — > fib via S R (Vq 1 ) for M LQ < y/s supersede previous 



upper bounds [0 from B -> /ie, for M LQ below 200 GeV (220 GeV). On the other 



hand, the limits from /i conversion in titanium [jj and from forbidden K decays 
which involve only first and second generation quarks are much stronger than the 
corresponding ZEUS limits. 



For M LQ below 200 GeV, the ZEUS limits on ed 



previous limits from r 
and b, respectively. 



ne 



35 



rg 2 through Sq and V R supersede 
for q 2 = d, s, 



Ke 36 , and B -> reX 37 



Figure |7| illustrates examples in which the existing low-energy limits, though less stringent 
than the ZEUS limits at low M LQ , become more stringent at higher masses. As described 
in the next section, this is not always the case. 

An alternative approach to setting limits, which was employed in reference [0 is to assume 
that the branching ratio B iq2 is givenQ by X 2 lq2 /(X 2 eqi + Xj q2 ), and to set limits on Xi q2 for 
a fixed value of X eqi . Such limits are shown in figure |8|. For F = LQs, our limits are 
similar to those of reference [0, while for F = 2 LQs, the ZEUS limits are stronger due 
to inclusion of e~p data. 

11 This formula for Bi q2 assumes that the leptoquark does not couple to neutrinos. 
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Finally, a third way to illustrate the sensitivity is to assume that the LQ couplings have 
electromagnetic strength (\ 2 eqi = \ 2 m2 = 47ra for LQs which couple to e and /i) and to 
determine a lower limit on the allowed LQ mass. Such limits are shown in table |l]. For 
scalar leptoquarks, lower mass limits between 207 GeV and 259 GeV are set. Somewhat 
stronger mass limits, between 219 GeV and 272 GeV, are set on vector leptoquarks for 
which both the production cross section and detection efficiency are higher. 



5.2 High-Mass Leptoquark Limits (M LQ > 300 GeV) 

For high-mass leptoquarks, the cross section is proportional to the square of ^ eq ^ q , 2 = 
\ eqi \e q2 /M 2 Q , and a factor which does not depend on either the leptoquark couplings 
or mass (see equations || and This is also true of rates for lower energy forbidden 
processes For a given limit on ^/ eqi e q2 , the limit on the product A egi A£ g2 is proportional 
to M\ Q . As M LQ increases, the upper limit on the product of the couplings eventually 
exceeds unity and the perturbation expansion, on which the formulae in the appendix are 
based, breaks down. Even so, the parameter ^ eqi eg 2 serves as a reasonable figure of merit 
for experimental comparisons. 

Tables 0, [5|, |] and ^| summarize the 95% CL upper bounds on ^ eqiiq2 , in units of 1CT 4 
GeV -2 from this experiment and also from previous experiments [10|. Here q\ and q2 



are the generation indices of the quarks which couple to e and to £ respectively^- Two 
important characteristics of these tables are summarized below. 



1. In the e \i case, for LQ species Vy 2 , Vy 2 , Sy 2 , or S^ 2 , the limits from this exper- 
iment supersede prior limits in some cases where heavy quark flavors are involved, 

2. For the e <-> r case, we also improve upon existing limits for the same LQ species 
as in point 1. In addition, because the existing limits on e «-> r are much weaker 
than those for e «-> //, the ZEUS limits are the most stringent for several additional 
LQs which couple to c or b quarks. 



5.3 Limits for ftp Squarks 

Coupling limits for So and S\/2 leptoquarks were converted to coupling limits on d, u, 
and t as described in section |2]^. Figure ^] shows 95% CL limits on coupling vs. mass for 
ftp squarks which decay to \iq and rq. Here we assume the couplings at the production 
vertex (A' llfc for d k , Ay X for £t J ) and at the decay vertex (A^- fc ) to be equal. The solid curves 
are the ZEUS limits which are given for two assumptions. The lower curves (and the t 
limits) assume that all squark decays are i?-parity violating. The upper curves illustrate 
the impact of gauge decays on the limits. They assume that a single i?-parity conserving 
decay, namely q — >• q^ exists and that the photino is much lighter than the squark. Since 
our analysis is not sensitive to such decays (for which the branching fraction is given by 
equation |6|) these limits are somewhat weaker. The stop mixing angle is assumed to be 

12 Certain entries in these tables have been corrected and/or updated from reference jlO) after consul- 
tation with the authors j38| . 
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cos 2 9t = 0.5. The dashed curves are limits from low-energy experiments, adapted from 
reference |fL0"| . Table |6| gives lower mass limits for d, u, and t assuming that the couplings 
at the production and decay vertices are equal to the electromagnetic coupling (y/Ana). 
As with the low-mass leptoquark case described earlier, the ZEUS limits improve on 
existing limits in cases where quark flavor change accompanies the lepton flavor change, 
especially for e <-> r flavor changes. 



6 Conclusions 



We have searched for signatures of lepton-flavor violation with the ZEUS detector. Hy- 
pothetical exotic particles such as leptoquarks could induce lepton-flavor violation ob- 
servable at HERA. The tight constraints from sensitive searches for processes such as 
muon conversion in titanium and rare muon and meson decays do not apply to all pos- 
sible cases of LFV, many of which could be seen in ep collisions. Using 3.8 pb~ l of data 
taken at HERA during the 1993 and 1994 running periods, we have found no candidate 
events for LFV. The data permit us to constrain specific leptoquark coupling strengths 
as small as 10~ 3 a and to exclude leptoquark masses as large as 270 GeV (for electro- 
magnetic coupling) with 95% confidence. For M LQ ^> y/s, we calculate upper limits on 
the product of lepton flavor violating couplings divided by the square of the leptoquark 
mass, \ eqi \e q2 / Ml Q , and directly compare these with existing bounds from rare decays. 
Especially for e <-> r flavor changes, ZEUS has improved on existing limits for many 
flavor-violating scenarios. 
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8 Appendix 



We summarize here the cross section formulae for processes involving leptoquarks which 
couple only to left-handed or right-handed leptons. Process (1) can be mediated by either 
s-channel or w-channel leptoquark exchange. For the s-channel process, eq\ — > £q 2 , the 
differential cross section, for unpolarized beams, can be written as: 



d2 ° = _J_ n (r a K Q M 2 s 2 x 2 f \ scalar LQ 

dxdy 32nxs qiy ' 1 (sx - M 2 LQ ) 2 + M 2 LQ T 2 LQ \ 2(1 - y) 2 vector LQ, y ' 1 

where q±(x, s) is the parton density^ for the initial state quark or antiquark, X egi and \i q2 
are the couplings at the production and decay vertices, and T LQ is the total width of the 
leptoquark. The partial width for decay into lepton £, and quark q 2 , is 



r> _ M \2 x | scalar L Q 



^ vector LQ, 



so that the typical LQ sought here has T LQ <C M LQ . In the narrow width approxima- 
tion, which holds when the variation of qi(x) is small as x is varied by 5x ~ T LQ /M LQ , 
integration of equation (|7|) leads to the formula ||, with B £q2 = Y iq2 /Y LQ . 

For the ^-channel process eq 2 — > iqi, the differential cross section is given by: 

^ - 1 q2(x , -„ ) r x { ^ - £ ( 9) 



etad?/ 327ras ' [ S x(l - ?/) + M 2 J 2 1 2 vector LQ. 



In the limit that M 2 Q ^> s, integration of equations (|7|) and @ lead to equations (§) and 
(|), which are accurate to better than 10% for M LQ > 500 GeV. 

Note that any leptoquark will take part in both s- and u- channel interactions. For 
example an Sq leptoquark will mediate the s-channel process e + u — > /i + c as well as the 
w-channel reaction e + c — > fi + u. 



For s- and it-channel processes, we have used s and — u respectively as the scale in the parton 
densities. If we had used Q 2 , the calculated cross sections would vary by less than 4% for initial-state u 
and d quarks and by less than 16% for initial-state s, c, or b quarks. 

on 
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Table 1: The 95% confidence level lower limits on the LQ mass (GeV) for nominal elec- 
tromagnetic coupling. For leptoquarks which couple to /i, we set \ 2 eqi = A^ 92 = 4ira. For 
leptoquarks which couple to r, we set \ 2 eqi = \ 2 q2 = Ana. Limits are shown for all scalar 
(S) and vector (V) leptoquark species. 
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Table 2: The best upper bounds on A egi A Mg2 /Mf Q for F = 2 leptoquarks, in units of 1CT 4 
GeV -2 . Each column corresponds to a given leptoquark species and each row to the quark 
flavors qi and q<i which couple to e and fi, the generation indices of which are specified in 
the first column. The top line in each box gives the previous measurement JTDJ which had 



obtained the strictest limit. The limit from that experiment is given on the second line in 
the box and the ZEUS limit, shown on the third line, is printed in boldface if it supersedes 
the previous limit. The asterisks denote those cases where lepton flavor violation occurs 
only via processes involving top. 
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Table 3: The best upper bounds on A egi A Mg2 /Mf Q for F = leptoquarks, in units of 1CT 4 
GeV -2 . Each column corresponds to a given leptoquark species and each row to the quark 
flavors qi and q<i which couple to e and fi, the generation indices of which are specified in 
the first column. The top line in each box gives the previous measurement JTT1 which had 



obtained the strictest limit. The limit from that experiment is given on the second line in 
the box and the ZEUS limit, shown on the third line, is printed in boldface if it supersedes 
the previous limit. The asterisks denote those cases where lepton flavor violation occurs 
only via processes involving top. 
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Table 4: The best upper bounds on X eqi \ TQ2 / M\ Q for F = 2 leptoquarks, in units of 1CT 4 
GeV -2 . Each column corresponds to a given leptoquark species and each row to the quark 
flavors qi and q% which couple to e and r, the generation indices of which are specified in 
the first column. The top line in each box gives the previous measurement [TD[ which had 



obtained the strictest limit. The limit from that experiment is given on the second line in 
the box and the ZEUS limit, shown on the third line, is printed in boldface if it supersedes 
the previous limit. The asterisks denote those cases where lepton flavor violation occurs 
only via processes involving top. 
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Table 5: The best upper bounds on X eqi \ TQ2 / M\ Q for F = leptoquarks, in units of 1CT 4 
GeV -2 . Each column corresponds to a given leptoquark species and each row to the quark 
flavors qi and q% which couple to e and r, the generation indices of which are specified in 
the first column. The top line in each box gives the previous measurement [TD[ which had 



obtained the strictest limit. The limit from that experiment is given on the second line in 
the box and the ZEUS limit, shown on the third line, is printed in boldface if it supersedes 
the previous limit. The asterisks denote those cases where lepton flavor violation occurs 
only via processes involving top. 
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Table 6: 95% CL mass limits (GeV) for squarks with /^-couplings, of electromagnetic 
strength (A? lfc = Xf jk = Ana = An/ 128). The first line gives the mass limits for the case 
where the total branching fraction for gauge decays is 50%, the second line gives the mass 
limits for a squark which always has fip decays. The limits which assume 50% gauge 
decays are weaker because we did not search for the gauge decays. The mixing angle of 
the stop is assumed to be cos 2 9 t = 0.5. 
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Figure 1: The (a) s-, (b) u—, and (c) t-channel Feynman diagrams for LFV. For the 
s-channel and it-channel diagrams, we denote the couplings as Xe q , where the indices refer 
to the lepton and quark flavors. 



on 




Figure 2: Existing 95% CL limits |1(J on the product of the couplings \ eqi \t q2 vs. lepto- 
quark mass M LQ for a Sq leptoquark mediating e <-> \i and e <-> r transitions accompanied 
by various quark flavor changes. Each limit curve excludes the region above it. The verti- 
cal lines indicate lower limits of allowed leptoquark mass from LEP jy] and the Tevatron 



T2|j . The Tevatron limits apply to any scalar leptoquark which couples to /i and depend 
on the branching fractions to the [iq final state. The limits are shown for equal to 
0.5 and 1. 




c) 

Figure 3: Rp violating single squark production in ep collisions. Diagrams a) and b) 
show production of u and d squarks with leptoquark-like Rp decays, where t denotes the 
final-state charged lepton of generation i. The indices j and k denote the generations of 
up-type and down-type (s) quarks respectively. Diagram c) shows u production with an 
i?p-conserving decay. 
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Figure 4: Net transverse momentum (/ > t ca/ ) distribution of events in the class ^ sample (a) 
and the class e sample (b). The points represent the data. The solid lines show the Monte 
Carlo prediction, which includes CC DIS, NC DIS, resolved and direct photoproduction 
and 77 interactions. In the top plot, the shaded region shows the Monte Carlo prediction 
for all processes except CC DIS. In the bottom plot, the shaded region shows the Monte 
Carlo prediction for all processes except NC DIS. 
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Figure 5: The upper limit on the coupling at the production vertex (X eqi ) times the square 
root of the branching fraction to the \iq or rq final state (Bi q2 ) vs. leptoquark mass M L q, 
at 95% CL for scalar leptoquarks. The horizontal line indicates nominal electromagnetic 
coupling (Ag = Ana = Atc/128) and B lq2 = 0.5. 
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Figure 6: The upper limit on the coupling at the production vertex (A egi ) times the 
square root of the branching fraction to the \iq or rq final state (Be q2 ) vs. leptoquark 
mass Mlq, at 95% CL for vector leptoquarks. The horizontal line indicates nominal 
electromagnetic coupling (Ag = Aira = An/128) and B tq2 = 0.5. 
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Figure 7: The 95% CL upper limits on X eqi vs. leptoquark mass Mlq, for selected LQ 
species which decay to £q, where £ = \i (above) or r (below), assuming B iq2 = 0.5. The 
solid curves are ZEUS results and the various broken curves show existing limits Jl(J . 
Paired numbers in parentheses indicate the generations of the quarks which couple to e 
and £ respectively. 
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Figure 8: The upper limit on the coupling at the decay vertex (\e q2 ) vs. leptoquark mass 
M L q, for several values of the first-generation coupling at the production vertex (A egi ). 
Each curve is labeled by the value of \ egi . The dotted curves are for F = 2 leptoquarks 
and the solid curves are for F = leptoquarks. 
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Figure 9: Limits on the ftp coupling A- jfc at 95% CL for squarks which decay to fiq 
(i = 2) or rq {i = 3). For d limits, we assume that A' llfc = A- jfc , while for u (including 
t) limits, we assume that A' Xjl = A- jfc . The lower solid curves give the ZEUS limits for 
squarks which decay purely via i?-parity violation. The upper solid curves give the ZEUS 
limits for the case where, in addition, the gauge decay q — > exists, and where the 
photino is much lighter than the squark. In the case of t, we only consider ftp decays 
and assume that the stop mixing angle is given by cos 2 6 t = 0.5. The values of j and k 
for which each limit applies are indicated above the plots. The dashed and dotted curves 
show the limits from other experiments (adapted from [[10]). These limits depend on the 
generation of the quark which couples to the [i or r as indicated. 



Q7 



